Disorder effects during the evolution from BCS to BEC superfluidity 



Li Han and C. A. R. Sa de Melo 
School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, 

(Dated: April 27, 2009) 



USA 



O 
O 

(N 

Ph. 
< 

(N 



C/2 



I 

o 



> 



o 

ON 

o 



X 



We describe the effects of disorder on the critical temperature of s-wave superfluids from the 
BCS to the BEC regime, with direct application to ultracold Fermi atoms. In the BCS regime 
the pair breaking and phase coherence temperature scales are essentially the same allowing strong 
correlations between the amplitude and phase of the order parameter. As non-pair breaking disorder 
is introduced the largely overlapping Cooper pairs conspire to maintain phase coherence such that the 
critical temperature remains essentially unchanged. However, in the BEC regime the pair breaking 
and phase coherence temperature scales are very different such that non-pair breaking disorder can 
affect dramatically phase coherence, and thus the critical temperature, without the requirement of 
breaking tightly-bound fermion pairs simultaneously. Finally, we find that the superfluid is more 
robust against weak disorder in the intermediate region between the two regimes. 



PACS numbers: 03.75.Ss, 03.75.Hh, 05.30.Fk 

Ultracold atoms are special systems for studying su- 
perfluid phases of fermions or bosons at very low tem- 
peratures, because of unprecedented tunability. In par- 
ticular, ultracold fermions with tunable interactions were 
used to study experimentally the so-called BCS-to-BEC 
evolution, and population imbalanced systems. In addi- 
tion, there are other interesting directions to be pursued, 
including studies of the BCS-to-BEC evolution in opti- 
cal lattices [l|, Q , and the effects of disorder during the 
BCS-to-BEC evolution, which would allow the very im- 
portant study of the simultaneous effects of interactions 
and disorder at zero [sj and finite temperatures 

In ordinary condensed matter (CM) systems the con- 
trol of interactions is not possible, and the control of 
disorder is very limited, because the disorder potential is 
not known and can not be changed at the turn of a knob. 
Thus, in standard CM the disorder is usually described in 
terms of defects or impurities, whose positions in the solid 
are assumed to be random. In ultracold atoms it is now 
possible to create controlled disorder using laser speckles 
or lasers with incommensurate wavelengths, which were 
used to study the phenomenon of Anderson localization 
in ultracold Bose atoms [1, 0] , but that could also be used 
to study disorder effects in ultracold fermions. 

Thus, here, we describe the finite temperature phase 
diagram of three dimensional (3D) s-wave Fermi super- 
fluids from the BCS to the BEC limit as a function of 
disorder, which is independent of the hyperflne states of 
the atoms and is created by a Gaussian-correlated laser 
speckle potential. Our main results are as follows. First, 
in the BCS limit the amplitude and phase of the order 
parameter are strongly coupled, such that pair breaking 
and loss of phase coherence occur simultaneously. In this 
case, the critical temperature is essentially unaffected by 
weak disorder, since the disorder potential is not pair- 
breaking and phase coherence is not easily destroyed in 
accordance with Anderson's theorem 7]. Second, in the 
BEC limit the breaking of local pairs and the loss of phase 



coherence occur at very different temperature scales. In 
this case, the critical temperature is strongly affected by 
weak disorder, since phase coherence is more easily de- 
stroyed without the need to break local pairs simultane- 
ously, and Anderson's theorem does not apply. Third, we 
flnd that superfluidity is more robust to disorder in the 
intermediate region between the BCS and BEC regimes. 

To investigate the physics described above, we start 
with the real space Hamiltonian {h = 1) density for three 
dimensional s-wave superfluids 

7i(x) = (x) f-Yl _ + Vd^si^)) 7A,(x)+J7(x), (1) 
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where J7(x) +/dx'T/(x,x')V'T(x')V'x(x')V'l(x)V'TW 
represents a term containing the interaction potential 
F(x, x') = —gS{x — x'), and V'J(x) represents the cre- 
ation of fermions with mass m and hyperflne state (spin) 
a. In addition, Vdis(x) represents the disorder potential, 
and /i, the chemical potential. We choose Vdis(x) to be 
independent of the hyperflne state, a choice that can be 
easily relaxed. 

Although there are many ways to introduce disorder, 
we will make a particular choice that the disorder po- 
tential is governed by a Gaussian distribution P[Vdis] = 
exp[— / dxV^^^{x.)/{2K)]/A, with normalization constant 
A = JV[Vdis]exp[- J dxy;f.^(x)/(2K)], which leads to the 
conflgurational average (Vdis(x)Vdis(x')) = k6{k — x'). 

To derive the effective action for a flxed conflguration 
of disorder, we deflne the local chemical potential /i(x) = 
/i— Vdis(x) and follow the functional integral formulation 
of the evolution from BCS to BEC superfluidity to 
obtain 5cff = 5o + ^g, where 5*0 is the action of unbound 
fermions in the presence of weak disorder and is given by 



5*0 [Fdi 
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/dx^ln [l 



e(k,x)/T 



(2) 



with ^(k, x) = Ck — /i(x), and ek = k"^ /2m. The sec- 
ond contribution to Sf.s corresponds to Gaussian pairing 
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fluctuations 

^G[Vdis] = ^J dx^ A(q)r-i(g,Fdis)A(<z), (3) 

where A(g) is the pairing field, and 

Xi + X2 



m^qi - Ci - 6 - 2Fdis(x))] 



C (4) 



is the pair correlation function in the presence of dis- 



order, where C 



-{mV/'iT:as) + Ek(2ek)" 



q 



(q, iqi) represents the four-momentum, the function 
Xi = tanh [(^1 -I- Vdis(x)) /2T] describes the occupation 
of fermions with energy ^1 — ^(k— q/2), and the function 
X2 = tanh[(^2 + ^is(x)) /2T] describes the occupation 
of fermions with energy ^2 = C(k + q/2). 

The description above corresponds to a semiclassi- 
cal approximation which is valid for weak disorder po- 
tentials. However, quenched (as opposed to annealed) 
disorder is characterized by a static disorder potential 
and the necessity to average the thermodynamic poten- 
tial ri(Vdis) rather than the partition function Q. The 
configurationally averaged thermodynamic potential is 
^avi^,) = {^{Vdis)), and depends on disorder via the 
parameter n. The thermodynamic potential for a fixed 
configuration is f2(Vdis) = —ThiZ(Vdis) when expressed 
in terms of the partition function Z{Vdis) = ^o(Kiis) x 
•^G(Vdis), where Zo(Kiis) = exp [-5*0(^118)] is the par- 
tition function for unbound fermions, and ZciVdis) — 
J T>[^, ^] exp [-S'g(^, Vdis)] is the partition function 
for the pairing field. Thus, il(Vdis) = f^o(Vdis) + f^G(Vdis), 
where the unbound fermion thermodynamic potential is 
f^o(Vdis) = -TlnZo(Vdis) TSo{Vdis) while the pairing 
field contribution is ricCKiis) = ^TlnZG{Vdis), which is 
approximated by 



^G{Vdis) = ~^ J dx^ln[Tr(g,ydis(x))] 



(5) 



Expanding in Vdis and taking the configurational av- 
erage leads to f2o,av(K) = T{So[Vdis]) = f^o,av(0) -I- 
AJ7o,av(K) for the thermodynamic potential of un- 
bound fermions, and ^Cavi^) = ^G,aviO) + ^^G,av{n) 
for the thermodynamic potential due to the pair- 
ing field. Here, no,av(0) and flcaviO) are just the 
thermodynamic potentials without disorder |8|], while 
the disorder-dependent corrections to the thermody- 
namic potential are Ar^o.avC's) = "^^^k L{\s.)k, with 
L(k) = 1/ [2T2(l + cosh(a/'r))], and AfiaM^) = 



-TJ2„Miq)K with M{q) 



P{q) - |g(g)|V2. 



The first term of M (q) can be expressed as 

Piq) = EkELo Cn [^r + ^2"] ^o{q)/V where Toiq) = 
r{q,Vdis = 0) and = tanh [(C™)/2r]. In addition, 
the coefficients found in the expansion are cq — — C2 = 
l/(2rz), ci = (16r2 - z^)/{8T^z^), and C3 = l/{8T^z), 
with z = iqt — ^1 — ^2- The second term is Q{q) — 



J2wJ2l=o'^n[Xl' + X^]To{q)/V where the coefficients 
are do = -^2 = l/(4Tz), di = 1/z^. 

The total thermodynamic potential is then f2av(K) = 
f^o,av('i) + f^G,av('«), and the corresponding number equa- 
tion is = — 9riav(K)/9/i, where 



N ^ No{K) + NGiK), 



(6) 



with Nq{k) = dflo,av{K)/dii being due to unbound 
fermions, and Ng{k) = dflG,av{K)/d^ being due to pair- 
ing. In the limit of k ^ 0, the standard results for the 
thermodynamic potential, number equation, and time- 
dependent Ginzburg-Landau theory are found 8]. How- 
ever, in the presence of disorder and in the vicinity 
of Tc, the low-energy and long- wavelength Lagrangian 
density Cg — J~-nd + derived from the action of 
Eqs. ([2]) and ([3]) has two contributions. The first is a 
non-dissipative part 



dr 



2rn, 



- fi^ -t-7Vdis(x) 



where ^ = 4'(x, r). Here, the term containing the effec- 
tive mass m* is the kinetic energy of the pairing field, 
fi^ plays the role of the pairing field chemical potential, 
is the effective interaction, and 7 Vdis (x) is the effec- 
tive disorder potential. The second term is the dissipative 
contribution reffecting the decay of fermion pairs into un- 
bound fermions, which takes the Caldeira-Leggett form 



Zn 



dr 



,|vI/(x,T)-«'(x,r')| = 



Since we are considering the case of weak disorder, 
these parameters can be easily related to original param- 
eters of the theory without disorder 0] . Using the nota- 
tion X = tanh(,^k/2r) and Y = sech^(^k/2T), the orig- 
inal coefficients are a{^,T) — + Ek ~ 5^7 ' 
corresponding to the order parameter equation in the ab- 
sence of disorder when a{iJ,,T) = 0, 



X 



Y 



XY kl 



and d{fi,T) — d^ + idi, with dji — J^k^/'^^k and 
di = [TTN{eF)y/J^/{8Ty/ep)^ determine the effective 



mass m*, while b{^,T) = ^j^. 
mine the effective interaction g. 



X 



and d deter- 



For instance, the effec- 
tive mass is m* = d^jn/c^ the effective chemical poten- 
tial /i* — —a/du, the effective interaction is — b/d\, 
the dimensionless amplitude of the disorder potential is 
7 —{da/d^)/df>, and A = dj/dR. A natural choice of 
dimensionless parameters arc l/^kpo-s) for interactions, 
7] — Knp/ejp for disorder, and T = T/ep for temperature, 
where kp is the Fermi momentum, np = fcl/STT^ is the 
fermion density and cp — kp/2m is the Fermi energy. 
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Notice that the condition 

fi,{f,l/{kFa,))=0 (7) 

corresponds precisely to the order parameter equation in 
the absence of disorder a(/z, T) = 0, indicating that this 
equation is not expHcitly affected by weak disorder, as 
required by Anderson's theorem 0]. However, the deter- 
mination of the critical temperature Tc{r]) and the chemi- 
cal potential /Xc(77) as a function of dimensionless disorder 
ry for a fixed scattering parameter l/(fci?as) requires the 
simultaneous solutions of the number Eq. ([6]), and the 
order parameter Eq. ([7]). 
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FIG. 1: The pairing temperature Tp (dashed line) and the 
critical temperature Tc (solid line) for uniform superfluidity 
in units of ep as a function of disorder parameter rj on a) 
the BCS side for l/kpas — —1.8 and on b) the BEG side 
for l/kpcis = 1.2. The dotted line in b) corresponds to the 
temperature where the superfluid density vanishes. 

In the BCS limit, j —^ 2 indicating that the pairing 
field feels an effective disorder potential twice as large 
as that felt by individual unbound fermions, while A — 
tt/ (2Tc) ^ 1 indicating that the dynamics of the pairing 
field is highly overdamped due to the decay of fermion 
pairs into unbound fermions. In the BCS limit without 
disorder the order parameter equation alone = (or 
a = 0) determines the critical temperature @]. As seen 
above this equation remains unchanged in the presence 
of weak disorder. Thus, the critical temperature is essen- 
tially unchanged, because the fermion chemical potential 
fj, is very large, positive, and remains pinned to cf, in the 
case of perfect particle-hole symmetry. This is a manifes- 
tation of Anderson's theorem, and changes in the critical 
temperature can occur only via the disorder dependent 
shift in the chemical potential, which leads to a change 
in the single particle density of states. 

The number equation in the BCS limit is approximated 
by iV « No{0) + ANoi-q). The first coefficient repre- 
sents the contribution in the absence of disorder iVo(O) = 
^[^(1 — X), while the second contribution involves the ef- 
fects of disorder ANo{r]) = vJ2kXY/4:{f)^ . In the case 
of perfect particle-hole symmetry AiVo(?7) = 0, and in- 
deed = ep- Thus, strictly in the BCS limit {l/kpas 
—oo), where the assumption of perfect particle-hole sym- 
metry is made, the chemical potential is pinned to the 
Fermi energy (^ = e^) and the critical temperature is 



unchanged Tc = (e'''/7r)8e~^ exp(— 7r/2A;F|as|). Since the 
presence of disorder does not alter the total density of 
fermions np — kp /Stt^ , relaxing the condition of perfect 
particlc-hole symmetry makes the chemical potential ad- 
just itself to n = ep — Afi{ri), with Afi{ri) > leading 
to a corresponding reduction in the critical temperature 
Te(77) = r,(0) - AT,(77), with Ar,(7y) > 0. In Fig. [Di, 
the behavior of the pairing temperature Tp(rf) and the 
critical temperature Tc{rj) are shown in the BCS regime 
as a function of rj for X/kpag = —1.8. These curves are 
obtained by solving the corresponding number and order 
parameter equations for varying rj. 

As the attraction is increased towards unitarity and 
/i = 0, the relative change in Tc for fixed disorder pa- 
rameter 77, expressed as ATc/Tc[Q) = [I -Tc{ri)/Tc{Q)], 
decreases as the attraction strength is increased, and 
reaches a minimum as indicated in Fig. O Beyond this 
minimum ATc/Tc(0) increases again, being the largest 
in the BEC regime. These results indicate that super- 
fluidity is more robust to disorder in the intermediate 
regimeof — 1 < l/(fci?Os) < 1, where the coherence length 
reaches a minimum Q (or the critical current reaches a 
maximum). However, the relative change ATp/Tp{Q) = 
[1 — Tp{r\)lTp(^)\ for fixed disorder always decreases with 
Xjikpag) indicating that it becomes increasingly more 
difficult to break pairs with larger binding energy as the 
BEC regime is approached. As discussed next, the dif- 
ference in behavior between ATclTc{^) and ATpjTp^G) 
in the BEC regime is attributed to temporal phase fluc- 
tuations, which become increasingly more important as 
Xjikpa^^ increases. 
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FIG. 2: The relative change Arc/rc(0) (solid line) and 
Arp/rp(0) (dashed line) throughout the BGS to BEG evolu- 
tion for r\ — 0.1. Notice the minimum in ATc/Tc{0) occuring 
in the intermediate region —1 < l/kpas < 1 indicating that 
the superfluid is more robust to disorder in this region, while 
the Arp/rp(0) always decreases indicating that the pairing 
temperature becomes less sensitive to disorder in the BEC 
regime. The dotted line corresponds to the analytical result 
in the BEG limit. 

In the BEC regime, 7-^2 indicating that the pairing 
field feels an effective disorder potential twice as large as 
that felt by individual unbound fermions, but A ^ indi- 
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eating that the dynamics of the pairing field is undamped, 
reflecting that tightly-bound fermion pairs have very long 
lifetimes. In the BEC regime {l/kpag ^ 1) the number 
equation is dominated by the contributions coming from 
Nq 3> Nq given that all fermions are paired into molec- 
ular bosons. Here, Nq is well approximated by the num- 
ber of bound fermions (molecular bosons) A'b, since the 
number of scattering states A^sc is very small. The con- 
figurational average can be performed using the replica 
technique 0, [13] on the Lagrangian density Cg lead- 
ing to the action Sg = {'^/T)Y.qb{q)G-^{q)h{q), with 

G^'^{q) = M* ^ ^ ^(9)1 where the self-energy is 

IkKlk.l ^ 

Here, G{q) plays the role of the molecular boson propa- 
gator and leads to the thermodynamic potential fJav ~ 
nG,ecv = ~2TY,q[^-a-G-'^{q)/T] with total number of 
fermions N w 2iVB, where Nb = TY^^C^^iq), and 
renormalized chemical potential /i* — fi^ + g^iTB + 
7^'*S|k|<|k I Upon summation over Mat- 

subara frequencies iqg, we arrive to the boson density 
ub = Nb/V corrected by disorder 

UB = C(3/2)(m,r,/27r)3/2 ^ ^T,ml / in^ , (8) 

when [1^ = 0. The solution for Tc can be ob- 
tained numerically, but an analytical solution is pos- 
sible when T] = kuf/^p ^ 1- This leads to 
Te = Tc(0) [1 - -f^KT,{0)ml/6Tr^nB] , where T,(0) = 

the Bose-Einstein condensation 
temperature for a gas of molecular bosons, which be- 
comes Tc{0) — 0.218ef^, when — + 2m and hb np/2. 

Unlike in the BCS regime, the critical temperature Tc 
is directly affected by the presence of disorder in the 
BEC regime, and is determined essentially by the num- 
ber equation. A small amount of disorder affects the 
phase coherence of the molecular bosons through the 
emergence of an incoherent part in the molecular bo- 
son Green's function G^^{q) manifested by the branch 
cut -^/l/l* I — iqi^ present in the self-energy £((7). If only 
the zero frequency {iqi = 0) contribution were con- 
sidered, one would be led to the conclusion that dis- 
order would not affect the critical temperature in the 
BEC regime. However, disorder introduces quantum 
(temporal) phase fluctuations in the BEC regime, where 
fermions are largely non-degenerate, particle-hole sym- 
metry is absent, and Anderson's theorem does not apply, 
thus leading to a strong suppression of Tc- 

The critical temperature calculated corresponds to the 
transition between the superfluid and a non-uniform 
phase containing local superfluid islands separated from 
each other by peaks of the disorder potential. The tran- 
sition from the local superfluid regime to the normal 
phase may be determined by imposing the condition that 



the superfluid density p(T, 77) — 0, which leads to the 
schematic phase diagram in Fig. [3l 




FIG. 3: A schematic finite temperature phase diagram as a 
function of dimensionless disorder rj and interaction 1/kpas 
parameters based on extrapolations for arbitrary disorder. 
Notice existence of regions of local superfluidity, without 
global superfluid order, which emerge at intermediate disorder 
before the local superfluidity is also destroyed by localization 
phenomena. Regions of Fermi- Anderson insulator (FAI) and 
Bose-Anderson insulator (BAI) are also indicated. 

We analyzed the effects of disorder on the critical tem- 
perature for superfluidity of ultracold fermions during 
the evolution from the BCS to the BEC regime. For s- 
wave superfluids, we showed that weak disorder does not 
affect the critical temperature of a BCS superfluid with 
perfect particle-hole symmetry in accordance with An- 
derson's theorem, as the breaking of fermion pairs and 
the loss of phase coherence must occur at the same tem- 
perature. However in the BEC regime, phase coherence 
is more easily destroyed by disorder without the need of 
simultaneously breaking fermion pairs. Finally, we also 
showed that the superfluid is more robust to the presence 
of weak disorder in the intermediate (crossover) region. 
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